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^ ■ Abstract 

\o '• 

, Coherent states provide an appealing method to reconstruct efficiently the 

' pure state of a quantum mechanical spin s. A Stern-Gerlach apparatus is used 

■ to measure (4s+l) expectations of projection operators on appropriate coherent 
Q>^ , states in the unknown state. These measurements are compatible with a finite 

I number of states which can be distinguished, in the generic case, by measuring 

r~| ' one more probability. In addition, the present technique shows that the zeroes 

Oh. of a Husimi distribution do have an operational meaning: they can be iden- 

I tified directly by measurements with a Stern-Gerlach apparatus. This result 

■ comes down to saying that it is possible to resolve experimentally structures in 
^ . quantum phase-space which are smaller than h. 

^ ^ ■ To determine the unknown state of a quantum system poses a problem which is inter- 

^ ■ esting from a theoretical, an experimental, and a conceptual point of view. To set up a 

^ . successful approach to state reconstruction (see and for reviews) means to iden- 

tify a collection of observables, a quorum, such that their expectation values contain 
all the information about the state. As for other inverse problems, "trial and error" 
is the typical approach since no systematic method is known to find an exhaustive 
quorum. The theoretician might focus on the search for a quorum easily realizable in 
the laboratory providing certainly a restriction and possibly an inspiration. On the 
experimental side, high standards in the preparation of individual quantum systems 
are needed since the measurement of expectation values requires reliable input states. 
By now, these standards have been met for quantum systems such as an electromag- 
netic wave 1^, vibrating molecules 0], ions caught in a trap 0, and atoms moving 
freely in space after scattering from a double slit [0. Conceptually, methods of state 
reconstruction offer exciting new perspectives since they allow one to faithfully repre- 
sent quantum mechanics in unconvential ways eliminating the wave function (or the 
statistical operator) in favour of positive probability distributions 0, |[ or expecta- 
tion values P]. In a similar vein, a quorum turns out to be an all-embracing source 
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of information allowing one even to determine the expectation values of operators for 
which no measuring apparatus is available |TU]. 



A brief summary of early work including the reconstruction of both pure and mixed 
spin states can be found in . Inspired by tomographic methods to reconstruct the 
density matrix of a particle, an entirely different approach has been worked out in 



12 , and it has been adapted in 113] to determine a single quantized cavity mode. 



Coherent states are used as a tool for state reconstruction of spin density matrices in 



14| , an approach which turns out to be relatively straightforward. 

Taking into account the extra information that the state to be reconstructed is 
pure — thus characterized by considerably less parameters than a mixed state — does 
not automatically simplify the problem. On the contrary, to eliminate the redundancy 
from the quorum makes it difficult to set up constructive methods of state determina- 
tion |jl5|. However, as shown in the following, coherent states offer a new and simple 
approach to determine a pure spin state by measurements with a Stern-Gerlach ap- 
paratus. 

The pure states of a spin of magnitude s belong to a Hilbert space Tig of 
complex dimension (2s + 1), which carries an irreducible representation of the group 
SU (2). The components of the spin operator S = hs satisfy the commutation relations 
[sx, Sj/] = "^Sz,..., and they generate rotations about the corresponding axes. The 
standard basis of the space Tig is given by the eigenvectors of the z component of the 
spin, Sz = Ta.z ■ s, and they are denoted by n^^), — s < yU < s. The phases of the 
states are fixed by the transformation under the anti-unitary time reversal operator 
T: Tl/i,!!^) = (— — /i, rix), and the ladder operators s± = Sx±iSy act as usual 
in this basis: 



s±|/i, n^) = ■y's(s + l)-/i(/i±l)|/i±l,n^). (1) 
With respect to the z basis one has 

l^)= E (2) 

fJ,= — S 

and the ray is seen to be determined by 4s real parameters since it is normal- 
ized to one and its overall phase is irrelevant. For the present purpose, however, a 
characterization different from (^) will be more convenient. 
Consider the eigenstates of the operator n ■ s, 

n ■ s l/i, n) = n) , —s<fi<s, (3) 

where the unit vector n = (sin "i? cosy), sin-i^siny), cos"!?), < < tt,0 < (f < 2tt, 
defines a direction in space. The probability Ps(n) to measure the value s with a 
Stern-Gerlach apparatus oriented along n equals 

Ps{n) = (s,n| V')(V'|s,n) = {^pl z ) {z\iIj) , (4) 

given thus by the expectation value of the operator \z){z\ projecting on a coherent 

\z) = |s, n) = exp[— i {} m((y9) ■ s ] |s, n^,) , (5) 



state !€. 
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where m{(f) = (— sin cos <y9, 0). In other words, the coherent state \z) is obtained 
from rotating the state \s, n^) about the axis m{ip) in the xy plane by an angle i). It 
is convenient to combine {'d^ip) into a single complex variable, z = tsin{'i)/2) exp[iip]. 
This provides a stereographic projection of the surface of the sphere to the complex 



plane. In terms of ^, a coherent state has the expansion |18 




, i.i.v. >: ,: z%s-kn,). (6) 



Project the state in (^) onto the coherent state 1^;): 

1 2s / r) \ 1/2 / * 2s 

providing thus a polynomial of degree 2s in z which has been factorized using its 2s 
complex roots z^, some of which may actually coincide. To specify their location in 
the complex plane requires 4s real parameters; therefore, they represent another way 
to characterize the pure state lip). It has been assumed in that the coefficient 
'?/'* ^ is different from zero which is true for almost all states. This "representation" 



of a spin state by 2s points in the complex plane [|T9[ has been used, for example, to 
study the quantum mechanical time evolution of spin systems possessing a classically 
non-integrable counterpart [^. Since the location of the zeroes z^ is easy to visualize. 



one has a convenient tool to study the quantum dynamics of the system through the 
paths traced out by the zeroes z^ in the course of time. 

The Husimi distribution 0], given by the squared modulus of (0), 
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is a non-negative function in the complex plane or, equivalently, on the classical phase 
space of the spin. Obviously, the Husimi distribution has the same zeroes as does the 
function (V'k). This observation allows one to implement a first method to reconstruct 
the spin state l^/'), a method which will turn out to be important from a conceptual 
point of view. According to Eq. (H), one just needs to finds those 2s orientations 
n„, n = 1, . . . , 2s, of the Stern-Gerlach apparatus for which the probability to detect 
the value s vanishes. Experimentally, this means to perform a two-parameter search 
which may be tedious but, in principle, is possible. Once the directions n„ and hence 
the zeroes z° have been determined, the state I?/') is known: multiply out the factors 
of the product in (|^) and you can read off the ratios V's-fc/V^-s- 

Before turning to a systematic approach to reconstruct the state it should 
be emphasized that, to the best of our knowledge, the operational meaning of the 
zeroes z^ exploited here has not been pointed out earlier. On the contrary, doubts 
have been cast on the physical meaning of the zeroes of the function P{z*, z) since one 
thought it impossible to access them experimentally, due to the coarse structure of the 
quantum phase space. The uncertainty principle for canonically conjugate variables 
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such as position and momentum was assumed to exclude the resolution of structures 
in phase space within an area of order %. For the spin, however, the vanishing of the 
Husimi distribution at a point in phase space can be tested in principle by means of 
a Stern-Gerlach apparatus. For a particle, a similar resoning is expected to apply. 
If one were able to measure the probability to find the system at hand in whatever 
coherent state, point-like structures in phase space such as the zeroes of its Husimi 
distribution would immediately be endowed with physical meaning. 

Here is a systematic approach to reconstruct a pure spin state based on similar 
ideas. It is a two-step procedure: first, the probabilities for the value s to occur are 
measured with a Stern-Gerlach apparatus oriented along (4s + 1) prescribed directions 
in space. These experimental data are compatible with 2^* pure states which can 
be written down explicitly. Second, by an additional set of mesurements along 2s 
directions, which are calculated on the basis of the first series of experiments, one can 
determine the spin state unambigously. Alternatively, the measurement of one single 
probability is sufficient in most cases. 

Step I: On the real G -R, the Husimi distribution turns into a real 

polynomial of degree 4s multiplied with a nonzero factor: 

\ib P "^^ 1 
^(^'^) = n 1 2^2. n - 2«n^ + = TT—^s E cao;" ' (9) 

(1 + x^Y' ^j-^ V ^ (1 + x^r' t^o 

where = Un + ivn, and the expansion coefficients are functions of the zeroes: Ca = 
C\{zi, Z2, . . ■ , Z2s)- Select now (4s + 1) distinct points on the positve real axis, 

Xi, = tan ^ > ' ^ Xiy' , < z/, z/' < 4s , (10) 

which correspond to (4s + 1) points n,y of the sphere on the great (semi-) circle with 
coordinates ipn = and 'dn following from (^). To fix ideas, one might pick 

u+l ^ 2s- u , , 

x^ = - — — cosi9^ = z/ = 0, l,...,4s, (11) 

4s + 1 — z/ 2s + 1 

providing directions with equidistant projections on the axis n^. Measure now the 
probabilities to obtain the value s along the directions rij^; this results in (4s + 1) 
numbers ^ ^ 

P{x,,x,) = p, = f^M,xCx, M,, = — ^p^, (12) 

A=0 U+^i/J 

with a real matrix M^x of dimension (4s + 1) x (4s + 1), the entries of which are all 
between and 1. Eq. (|^) can be solved for the coefficients c = (cq, Ci, . . . , 04^)"^ 
according to 

c=M-^p, (13) 



since the matrix M, being of type Vandermonde has a nonzero determinant, 

4, , ^ 1 xo ■■■ (xo)'^' 
detM =(117^ 



1 X4s ■■■ {X4s) 



4s 
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4s 



\u=Q V ' i^J / 0<u<u'<4:S 



since ^ Xyi ^ v ^ i/', according to (|TU|). An explicit expression of the inverse of a 
Vandermonde matrix is given in p3| . 

Knowing the coefficients c\ obtained with (|T3p, the function P{x,x) is determined 
through the sum in Eq. (j^), and one can work out its 2s factors (x^ — 2unX + 
This, however, is not sufficient to determine the state {ip): generically, there are 2^'^ 
distinct states giving rise to the same expansion coefficients c. To see this consider 
the state j G (1, . . . , 2s) with 
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w^ok) = 7f^(^-(^?r) n i^-^i), (15) 



r-sU) 

having thus the same zeroes as except for Zj which has been replaced by its complex 
conjugate (zj)*. The Husimi distributions P{z*,z) and Pj{z*,z) = \{ip{j)\z)\'^ differ 
from each other for general z but for real arguments (cf. Eq. (^) they coincide: 

P{x,x) = P,{x,x), jG(l,...,2s), (16) 

since the information about the sign of the imaginary part of the j-th zero has dropped 
out. This is also obvious from looking more closely at the factors [x"^ — 2unX + I^^P) 
in Eq. (|^): they allow one to calculate unambiguously the real parts m„ of the zeroes 
but only the modulus of their imaginary parts, = — |u„p. Thus, a twofold 
ambiguity is present for each of the 2s zeroes of\tp), giving rise to a total of 2^* statesQ 
compatible with the data p. It is straightforward to remove the ambiguity by a second 
series of measurements. 

Step II: It remains to find out whether a zero is located at z^ or at (z^)*, n = 
1, . . . , 2s. To do this, one orients the Stern-Gerlach apparatus along the axis n„ or 
along n„. associated with the point 2° or its complex conjugate, respectively. If the 
probability to measure the value s along n„ does not vanish, the zero in question is 
located at n„*. If it vanishes, the zero is determined by the actual orientation of the 
apparatus. Repeating this procedure for the 2s pairs of points, one is able to remove 
the ambiguity completely. 

There is a possibility — working in almost all cases — to reduce the measurements 
required for Step II to the determination of a single probability. The 2^* states compat- 
ible with the data p are known explicitly. Now calculate the values of the probability 
to measure the value s along one randomly chosen direction n' in space which does 
not coincide with any of the putative zeroes. If the 2^* calculated values all differ from 
each other, it is sufficient to perform the measurement along n' in order to identify 
the correct state. If, however, two or more probabilities with respect to the axis n' 
coincide, a slight change n' n" usually will be sufficient to remove the degeneracy. 
Unfortunately, it is difficult to actually determine a particular spatial direction such 



""^This ambiguity happens to be smaller in exceptional cases: one or more of the zeroes might be 
real, or some zeroes might coincide. 
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that the probabihties are guaranteed to have different values — although this is the 
generic case to happen. 

When carrying out the first proposal one measures a total of (4s+ 1) + 2s = (6s + 1) 
real parameters in order to determine the state {ip). The second method requires only 
(4s + 1) + 1 = 2(2s + 1) measurements which exceeds the number of free parameters 
by two only. It is essential, however, that both cases give results linear in s providing 
thus a substantial reduction compared to the optimal method (oc s^) for a mixed 
state. 

Points Xj, on the positive real axis are not the only choice for a successful recon- 
struction scheme. They might be located, for example, on a half line starting at the 
origin which encloses an angle G (0, 27r) with the real axis. This situation is ob- 
tained from the original one by a conformal transformation, z ^ w = zexp{—iip), 
that is, a rigid rotation about the origin by an angle if. 

It is also possible to carry out Step I with (4s + 1) distinct complex numbers of 
modulus one, exp[i(y9j,], providing directions n^, on the equator of the sphere. The 
elements of the matrix M then read 

M., = i^(exp[z^,])' , (17) 

its determinant being different from zero if no two angles ipi, coincide. Again, a 2^*- 
fold ambiguity remains which is due to the existence of other states with the same 
Husimi function now on the unit circle. Consider, for example, the state 

Its Husimi function for the argument z = exp[iif] satisfies 

P,(e-^^ = Pie-'"', e'^ , j G (1, . . . , 2s) , (19) 

so that one cannot distinguish the states {ipij)) and {ip) by looking at the unit circle 
only. The property (p!9| ) is due to the fact that P{z*, z) on the unit circle is invariant 
under the transformation z — > 1/z*, while the invariance associated with ( [T6|) is due 
to a reflection z — > z* about the real axis. 

Further, one can displace the points rigidly, z—>-w = z + z,za. fixed complex 
number, subsequently reworking one's way through the derivation given above. A shift 
of the points on the positive real axis to the left shows once more that one can use 
also (4s + 1) points on a (and therefore: any) great circle. Apply now a shift z to the 
ensemble of points located on the equator: then, this great circle turns into another 
circle on the surface of the sphere providing thus another ensemble of directions for 
the measurements. 

It appears plausible — although no proof is available yet — that Step I can be based 
on a set of (4s + 1) arbitrary distinct points in the complex plane. One might speculate 
that such a generalisation is more effective than the method presented above: possibly 
no ambiguity will show up if (4s -|- 1) 'random' points in the complex plane are 
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chosen. In other words, the necessity of Step II (to perform additional measurements) 
could result from the non-generic choice of the points on a line (or circle) introducing 
unintentionally a reflection (or inversion) symmetry. 

Summing up, the method of state reconstruction via coherent-state expectation 
values can be based on any (4s + l) points in the complex plane located on an arbitrary 
circle, including its limit of a straight line. In the generic case, one additional mea- 
surement is needed to determine the actual pure state of the system. Furthermore, for 
a pure spin state, a search for the zeroes of its Husimi distribution provides another 
strategy of state reconstruction since they are seen to be accessible in experiments 
with a Stern-Gerlach apparatus. 
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